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ABSTRACT 

This paper is devoted to Radial Orbit Instability in the context of self-gravitating 
dynamical systems. We present this instability in the new frame of Dissipation- Induced 
Instability theory. This allows us to obtain a rather simple proof based on energetics 
arguments and to clarify the associated physical mechanism. 
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1 INTRODUCTION 

Instabilities in self-gravitating systems are fundamental pro- 
cesses to understand the shape and physical properties of 
objects such as galaxies or globular clusters. So far, only a 
few of such mechanisms are described in literature, namely 
Jeans instability, which governs the collapse of homogeneous 
systems; gravothermal catastrophe, which concerns isother- 
mal spheres; and radial orbit instability, which occurs in 
anisotropic, strongly radial spherical systems. If the first 
two are well understood, and have tak en their place in the 
study of dynamical stellar systems (see lSinnev fc Tremaind 
Il987 . sections 5.2 and 7.3), as a fact, the situation of ra- 
dial orbit instability is less clear. A complete story of this 
physical process, spanni ng almost forty years, is presented in 



pnysical process, spanni ng almost torty years, is presented m 
iMarechal fc Pere j jiooa V Three main points stick out (see 
the review for all detailed references): there is as yet no sim- 
ple analytical proof of this phenomenon; there is no global 
consensus about its actual physical mechanism; and yet it is 
a fundamental process which affects the phase space distri- 
bution of primordial gala:xies and contributes to produce the 
radial density profile of evolved systems. The present paper 
will address the first two of these points. 



1.1 Collisionless Boltzmann— Poisson System 

We consider a system constituted of a large number of 
gravitating particles interacting together. We will assume 
that all those particles have the same mass m. We denote 
as q and p the position and the associated impulsion of a 
particle with respect to some Galilean frame TZ, and T = 
(q, p) the corresponding point in the phase space R^'. 

We assume that the statistical state of the system is de- 
scribed at each instant t by a distribution function / {T,t), 
with / (r, t) dr representing the number of particles con- 
tained in the elementary phase space volume dF located 
around F. 



If the influence of collisions on the overall dynamics is 
neglectec(3, this distribution function solves the Collisionless 
Boltzmann-Poisson system (hereafter CBP) 
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at 



+ ^ ■ Vq/ - mVqV ■ Vp/ = f + {/, E} = 
VqV = 47rGm / /dp 



with boundary conditions ip =|q|^+oc O (r ^) and 
lim f — 0. The function tp (q, t) is the gravitational 

|q|,|p|-» + oo 

potential created by the particles, 



S(q,p,t) 



2m 



-f- mip 



is the one-particle Hamiltonian, and {., .} denotes the Pois- 
son Bracket defined by 

{/l, M = Vq/l • Vp/2 - Vq/2 ■ Vp/l 

Any stationary solution fo (F) of the CBP system is as- 
sociated to an equilibrium state of the particles distribution. 
It is now well known that CBP system is Hamiltonian with 
respect to a Poisson bracket of non-canonical form arising 
from the fact that a distribution function does not consti- 
tute a set of cano nical field variables (see iKandrupl Il990l : 
I Perez fc "Ahlll996l ): the set of distribution functions is an 
infinite-dimensional space. The total energy associated to a 
distribution function / can be written as 



H 



:/]=/dF|-/(F,.)4 



dF / dF 



,f{r,t)f{T',t) 



For any two functionals A [f] and B [f] of the distri- 
bution function, let {A, B) denote the Morrison bracket - 
introduced in the context of plasma physics by iMorris^ 



^ For a self-gravitating system with large values of A', this hy- 
pothesis is justified; see lBinnev fc Tremainel l ll987l) , part 1.2.1. 



2 L. Marechal, J. Perez 



l|l980l ) - defined by 



{A,B} = 



, ,SA SB] ^ 



where |j stands for the functional derivative of A, which is 
the linear part of ^4 [/ + 5f] — A [f] . One can easily obtain 
the Hamiltonian formulation of GBP system 

dF 



dt 



= {F,H) 



where F [f] is any functional of /. 



1.2 The stability problem 

The stability of equilibrium states is a very old problem of 
theoretical stellar dynamics, and a large variety of meth- 
ods has been used to tackle it. A clear consensus was found 
about the global stability of isotropic spherical systems with 
distribution function /p (E) monotonic ally decreasing: after 
the pioneering works bv lAntonovl(|l96ll '). linear stability was 
ob tained usin g energy methods after a long series of papers 
bv lKandrup fc Sv enet (1985) (and references within) or see 
also Perez fc Alv (jT99Q l for a comparison of the different 
results; using direct norma l mode techniques, complicated 
proofs were also obtained (|Fridmann fc Polvache"nkolll984l : 
|Palmeij[l994l ). Non- linear stability of such spherical isot ropic 
syste ms was also proven for some specific models f see iReinI 
|2002| . and references within). 

The stability of anisotropic spherical systems is a more 
difficult problem. The distribution function depends both 
on the one-particle energy E and on the squared one par- 
ticle angular momentum := p^q^ — (p • q )^. The most 
general result in this context was obtained bv lPerez fc Alvl 
l|l996t) and concerns linear stability for the restricted case of 
preserving perturbations, which includes radially symmet- 
ric ones. Non-linear stability is assured for some classes of 
generalized polytropes for whic h /p (E, L^) — E'^L^^ with 
adapted values of k and p (see iRein (|2002l ) and references 
within). 

Some very technical approaches using normal modes 
claim linear instability for anisotropic systems composed 
only of radi al orbits (see iFridmann fc Polvachenka 1 19841 : 
|Palmeijll99i): this is known as t he radial orbit instability. 

See lMerritt fc Agui lar' (^1985') for one of the first relevant 
numerical app roaches, Perez et al., ( 1996 ) for an intermedi- 
ate position or lBarnes et al.l ( 2009h and references within for 
the most recent situation of this problem. 

A co mplete historical account of radial orbit instability 
is given in lMarechal fc Perej (|2009 l). In the next section, we 
present some key features of this process. 



1.3 Basics of radial orbit instability 

Radial orbit instability (hereafter ROI) appears in self- 
gr avitating s ystem dyna mics wi t h the pionnering works 
of lAntonovl (Il973l) and iHenonI (|l973l ). A decade later 
iPolvachenko fc ShukhmanI (ll98lE propose a stability crite- 
rion based on the ratio of radial over tangential kinetic en- 
ergies, when it is to small the system must leave its spher- 
ical symmetry and fo r m a bar. This work is criticized by 
iPalmer fc Papaloizoul l|l987l ) which suggest, using normal 
modes techniques, that ROI can occur for arbitrary small 



values of the Russian ratio, provided the distribution func- 
tion of the system is unbounded for orbits with zero angular 
momentum. This paper is also the first one to propose a 
relevant physical mechanism to understand ROI based on 
resonant trapped orbits; we note that this mechanism needs 
a coupling between orbits. 

Several factors show that a radial system needs a "seed" 
from which R OI can appear. This is developed in detail in 
iRov fc Perej ([2004). The general idea is that there has to 
be a near-equilibrium state, so that coupling between or- 
bits has the time to develop, and the instability to grow. In 
iMac Millan et al.l ((1999) , density profiles in a power law are 
considered, which means that the core has the time to sta- 
bilize before outer zones collapse, causing ROI to appear; 
it also shows that adding clumps tends t o accelerate the 
process. A counterexample can be found in lTrenti fc BertinI 
(2006), which shows that homogeneous spherical haloes do 
not undergo ROI, as the system tends to isotropy before 
reaching equilibrium. For this reason, our study will focus 
on ROI emerging from equilibrium states. 

Although ROI is a natural candidate to produce tri- 
axiality w hich can occ ur in self-gravitating systems, it was 
noted (e.(7. lKatz|[l99ll ) that this spatial counterpart of ROI 
could disapear duri ng the merg i ng pr ocess of the galaxy 
formation. However, Huss et al.l (|l999l ) and more recently 
iMac Millan et all (|l999l ') have shown that ROI is a funda- 
mental initial process which shapes the phase space of the 
galaxy progenitor and allows it to get the good final mass 
density profile. It is therefore important to understand fully 
the nature of ROI. 

In this context, the objective of this paper is twofold. On 
the first hand, in section [2] we present a general method for 
investigate instability of self-gravitating syste ms. This ap- 
proach couples a general mathematical result bv lBloch et al.l 
( 1994) which generalizes Lyapunov theory, and the symplec- 
tic approach of the stability pro blem of CBP system (s ee 
iBartolomewlll97ll : lKandruDlll99(]|: andlPerez fc A.lvlll996l ) - 
it must be noted that KandruJ i 199 ll ) has already used this 
technique for non spherical systems without the complete 
mathematical background. 

On the second hand, in section |3] we apply this method 
to obtain a direct energy proof of the radial orbit instability 
when the system can dissipate energy. 



2 DISSIPATION-INDUCED INSTABILITIES 
AND SELF GRAVITATING SYSTEMS 

2.1 The method of energy variation 

Consider the first-order vari ation of an equilibr i um /p 
JV+/'^'- It is well-known (see \ 



tion ot an equilibr i um /p — » 
Bartolomewl (| 19711 ) . iKandrud 



1990) and iPerez fc Alvl (1199 6*)) that there exists a phase 
space function g{T,t), such that 

/W(r,o = -{3,/o} (1) 

This function is called a generator0 of the perturbation. 
Written in this form, j'^' is the largest class of physical 
perturbations which can be considered as acting on /p. In 
other words, f'^^^ is a deformation of /p and then there exists 



This function is clearly not unique. 
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a g such that we have equation ([T]). Associated to this per- 
turbation, variation of the total energy - which turns out to 
be of second order in q, see for instance a short calculation 
in lMarechal fc Pered l|2009i ) - is given by 



with 



{g,E}{g,fo}dr~Gm^ 



{.g,/o}{g',/o} 
|q- q' 



(2) 

where V refers to (q', p'), /g to fo (F') and so on. 

When _ff'^' [fo] is positive for a given set G of acceptable 
generators g, the system is reputed stable against the asso- 
ciated perturbations. This argument was detailed and used 
to prove, in the case when fo = fo (E) and dsfo ■= < 0, 
stability against all acceptable g; and when fo = fo [E,L^^ 
and dsfo < 0, stability for all g such that |q,L^| = which 
are called preserving perturbations (see IPerez fc Alvl 1 19961 
for all details). 

When there are negative energy modes, generators g 
that cause //'■'^^[/o] < 0, they are not necessarily associated 
to an instability. Taking into account dissipation in the sys- 
tem can drastically change its dynamics. 

In the next section, we will illustrate this point with a 
simpler, yet instructive example. 



2.2 An electromagnetic example of 
dissipation-induced instability 

Consider a particle of mass m — 1, charge e, let us de- 
note as q = {x, y, z)^ its position with respect to some 
Galilean frame. This particle is influenced by two forces: 
one derives from a potential V that is maximal at q = 0. 
We'll write V{q) = -^uj^q^. The other one is the Lorentz 
force generated by a static magnetic field B — Boe^ = VA A 
with A = {xBy — yex). The Lagrangian of this particle is 
JZ — ^q^ -f eq ■ A + ^uj^q^, then the impulsion p conjugate 
to the position q is given by p = Vq (£) = {px,Py,Pz)^ , 
and thus, with 13 = 



( px=±- Py 

< Py = y + Px 

y Pz= Z 

The Hamiltonian of the system is "H = p-q— £, the equations 
of motion are given by Hamilton's ones, i.e. q = Vp {%) and 
p = — Vq {H). The behaviour of {z,pz) being trivial and in- 
dependent from movement on the other axes, let us focus on 
the system in the reduced phase space of ^ = {x,y,Px,Py)^ 
for which one has 



5 = where A = 
with K 



ah PK 

1 

-1 



(3) 

(4) 



Since = — /2 one can see that if we split A — A + B 



with 



A = 




al2 



h 




and B ^ P 



K 
K 



we have the fundamental property AB — BA, so exp (At) = 
exp (At) exp (Bt) , by direct series summation one can find 
that 



exp (Bt) 



E(t) 



cos (pt) 
- sin (pt) 



sin iPt) 
cos (pt) 



G SO2 (I 



and 



drdr' 



exp (At) = 



aip-z (t) h 



V52 (t) h 



with 



VI it) = 

^2 (i) = 

•Pi (t) = 

f2 (t) = 



cosh (^at) 



if 



Q > 



Q < 



(-Q)"'''^sin (V-at^ 

The general solution of the problem then writes 

^ (t) = exp (At) ■ exp (Bt) ■ ^ (i = 0) 

and it is stable provided that a — uj^ — p'^ < 0. We note that 
this stability is not asymptotic as all eigenvalues of A lie on 
the imaginary axis. The mathematical condition on a corre- 
sponds to the physical case when the effect of the magnetic 
field B is stronger than the effect of the scalar potential V. 
We can see on this example that it is possible to have a sta- 
ble equilibrium even on a point where the potential is at a 
maximum; negative energy variations around an equilibrium 
IS not a sufficient criterion for an instability. The physical 
explanation is that the magnetic force, which does not de- 
rive from a scalar potential, tends to 'curve' the particle's 
trajectory, and if the magnetic field is strong enough, this 
can be enough to keep the particle close to the potential 
maximum in spite of the repulsive force. 

However, this behaviour is only possible as long as there 
is no energy dissipation. If the system is able to dissipate 
energy, such an equilibrium becomes unstable. For example, 
assume there is some form of fluid friction force F/ = — 7q, 
the system is no longer Hamiltonian. Movement is still trivial 
in the {z,pz) plane; keeping the same variables that we have 
used in the non-dissipative case, the equation of motion is 
now 



^ = Aj£, where A^, = A -I- ■yC, C = 





~PK 





-I2 



The matrix C happens to commute with B, so we have 
{A + 7(7)5 — B{A + 7C): the fundamental matrix of the 
dissipative system splits into 

exp (A-yt) = exp {[A + jC] t) ■ exp {Bt) 

As exp [Bt) is a rotation matrix, the stability of the dy- 
namics is governed by exp {[A -\- 'yC]t). The characteristic 
polynomial of A-y = A + yC is 

X(A) = A* + 27A3 + (7^ 



roots of which are 



A1.2 ~ 



and 



A3,4 = i [-7 ± Vl^ + 4a - 4i/?7j 



Let us focus on the transition from the stable equilibrium 
we have determined towards the dissipative case. We then 
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have a = cj^ — < and < 7 ^ 1. In this Umit case, one 
can get 

,2\ -1/2N 



and 

A3, 4 



-1± 1 



-1/2 



2\l/2 / ^ 

+0(7) 



From our assumption that a < 0, we have (1 — ^) > 1, 
therefore there is a pair of roots (eigenvalues of ^4^) with 
positive real parts. From the equation of motion it fol- 
lows that the system is unstable; this kind of insta- 
bility, linked to its operator's spectrum, is called a 
spectral instability. When 7 is not infinitesimal, this insta- 
bility persists as one can check by direct spectrum calcula- 
tion or by more elegant approaches. The physical meaning is 
clear: if the particle loses energy, the magnetic field cannot 
'curve' it back as close to the maximum as it was previously, 
and it will spiral further and further from the origin. 

2.3 Dissipation-induced instabilities 

The previous three-dimensional example is a special case 
of a general theorem which applies for finite dimensional 
systems: a Hamiltonian dynamical system with a negative 
energy mode (which could be stable without further hypoth- 
esis) becomes spectrally and hence linearly and non-linearly 
unstable when any kind of dissipation is introduced. This 
counterintuitive result takes its genesis from the classical 
works by Thomson (Lord Kelvin) and Tait (1879), but it 
was prov en only recently in the case of finite-dimensional 
syste ms (|Bloch et al.l |l994| and iKrechetnikov fc Marsdeiil 
I2OO7I ). and, as suggested by references in the latter, ap- 
pe ars to be very useful in mecha nics. More recent works 
bv lKrechetnikov fc MarsdenI (|2009l ') suggest that the infinite- 
dimensional case works similarly, although there is no defini- 
tive proof for the time being. In the context of theoretical 
astrophysics, it is interesting to note that H. Kandrup used 
such kind of arguments to inv estigate gravita tional insta- 
bilities for triaxial systems fsee iKandrupl flOQll) . before any 
actual, formal result. 

As recalled in section 11.11 CBP is a Hamiltonian 
infinite-dimensional system, so we can apply this theory of 
dissipation-induced instability for stability investigations in 
this context of gravitational plasmas. In the next section 
we will show that, when a spherical and anisotropic self- 
gravitating system becomes more and more radial, we can 
choose a certain class of g for which [fo] < 0: this proves 
the existence of negative energy modes in such systems. Fol- 
lowing the dissipation-induced instability theory such kind 
of gravitating systems will become unstable as soon as any 
kind of dissipation can appear. As noticed by Kandrup in his 
visionary paper, in physical self-gravitating systems dissipa- 
tion could take several forms like a little bit of gas, dynam- 
ical friction or at minimum gravitational radiation! In the 
context of numerical modelizations of self-gravitating sys- 
tems where radial orbit instability also appears, dissipation 
is also inevitably introduced by numerical algorithms of time 
integration or by potential computation. 



3 APPLICATION TO RADIAL ORBIT 
INSTABILITY 

3.1 Approaching a radial system 

A pure radial orbit system is characterized by particles with 
= 0, the corresponding distribution function could then 
be written /g" (i5, i^) — ip (E) S (L^) where is any pos- 
itive smooth normalized function, and 5 denotes the Dirac 
distribution. However, this distribution is very irregular in 
zero which is quite problematic, in addition to being unre- 
alistic (orbits can hardly be perfectly radial). So, instead of 
actually using the Dirac distribution, we will use functions 
that approach it. 

The choice we made is to use Gaussian functions. More 
specifically, we will consider an initial distribution function 
of the form 

fS {E, L') = ^ [E) [L') , 5.{L') = -L exp {-^ 

(5) 

By direct calculation one can easily check that, for any 
smooth function Z defined on the phase space, one has, by 
limited development of Z with respect to pe and p Q 



ZS.iL'W^ I ZS^L'^'^P^'^P' 



sin {9) 



r-2 1^+ 4 



-1-0(0 



Apr d q 



(6) 



which has a clear limit when a — > and selects the value of 
Z at L'^ = as expected. We could say that /q tends to a 
distribution function of purely radial orbits when a — > 0. 

In the following part, we will consider what happens for 
arbitrarily small values of a. 



3.2 Energy variation 

Our goal in this section is to show that there exist pertur- 
bation generators g that, for sufficiently small values of a 
(that is, for systems that are close enough to the purely ra- 
dial case), give a negative energy variation. To do so, we 
will start with a general g, calculate the energy variation 
[/q ] for our quasi-radial systems, and explain along the 
way what hypotheses we make about g to reach this goal. 



^ N.B.: variables pr, pg and thereafter are the conjugate vari- 
ables oir, 9 and 0, not the projections of p along the base vectors. 
We have pr = mr, pg = mr^9 and p^ = mr^ sin^{9)<j>. Also 



sm^{9) 

The calculation involves the well-known Gaussian integrals 



/ 



■ dx 



X e dx 



_ ^ 3 r 

c'^e dx = —r'^'jTv 
4 
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Usual Poisson bracket properties give, for fS{E,L^) 

{g, fS} = dEfS{9, E} + dL2fS{g, l^} 

where dsfS '■= and dL'ifS hence the second 

order energy variation ([2]) sphts into 



The second-order energy variation of perturbed near ra- 
dial orbit systems is then 



where 
Kl2 := 

Ke ■- 

SPL2 ■■ = 
SpE 



II 



{SpL2 +5pE){Sp',^ 



q'l 



-Sp'e) ,3 ,3 / 

d qd q 



G 



89 

SpL2Sp'j^2 



dpr d q 



L2=0 

dq dq' 



(21) 



(7) 

(8) 

(9) 

(10) 

(11) 



For a general perturbation, it is difficult to say more 
about the sign of H^^\ However, the system could receive 
any kind of perturbations. In order to go further, we have 
to make some assumptions about g. We already know, from 
section 11.21 that a radial function will not lead to an in- 
stability, and from section 12.11 that dependency on E and 
plays no part. To find a g function that works, we thus 
have to consider a non-radial perturbation. We can consider 
a perturbation that is axisymmetric around the z axis: 



dL2fS{g,E}{g,L^}dr 
dEfS{9,Efdr 

dL-2.fS{9,L^}<^'p 
m I dEfS{9,E}d''p 



g{r)=g{E,L\e,pg) 
With this hypothesis: 

{9,L'} 



2p,|| + 2p^^°^(^) ^3 



de 



He) dpe 



:{9,n 



With/,? = <^(£;)54L2^,weget 
dEfS = ^'{E)54L^) 



dL2fS 



(12) 

(13) 
(14) 

(15) 
(16) 



With the previous results, and using (|6}, we can calculate 
explicitly the four terms K]^2, Ke, 5pl2 and <5p_B in a power 
series of a for a — >■ 0. If we consider only the first term in 
o, which corresponds to the term of lower power in pe and 
P0, a long but straightforward calculation eventually leads 
to the following results: 



Kl2 = 



dpr d q 



L2=0 
2 



dpr d q 



SPL2 
5pE 



d'g 



+ 



L2=0 

cos{9) dg 



dOdpe sin(6) dpe 



dpr 



(17) 
(18) 
(19) 



L2=0 



ma I / , 



d^g cos{e) dg 



■ + 



dOdpg sa\{9) dpe 



dpr (20) 



L2=0 



As one can see Ke and 5pE are of order for a — >■ 0, 
whereas Ki^2 and Spi^2 do not depend on a in the same 
regime. Therefore, for near radial orbit systems one can ne- 
glect Ke in front of Kj^2 and SpE in front of 5pE2. 



The first term is clearly positive as an integral of a pos- 
itive function, while the second is clearly negative owing to 
the negativeness of the Laplacian operator: introducing 



M(q) ■- 

one has A/i = 4nSpE2 , hence 



Spl2&p'l2 ,a ,3 / 
^d qd q = 



Sp'l 



|q-q' 



■d q 



(22) 



■q'l 



_i_ 

47r 



/lA/idq 



The sign of i/'^' is thus unclear, unless we make another 
assumption about g. To be able to say more, we can consider 
a generating function verifying 



dg 



= while 



/ (23) 



For such perturbations one can easily check that Kj^2 
and 5pL2 7^ 0, therefore 



-G 



WK2^3q^3q,^Q (24) 



|q-q 

To summarize, in this section, we have obtained: 

Result 1. Let /o be a distribution function of nearly-radial 
orbits, such that 



fo = 'P (E) exp ( 



Let g be a generator of a perturbation, of the form (/(F) = 



g{E, L\ 9, pe), with ffj 



„ = while 



l2=o 



/O. 



Then for sufficiently small values of a, the energy vari- 
ation H^^^'I/q] caused by g is negative. 



3.3 Density variation 

It is interesting to analyse the density variation associated 
to the generator described in the previous section. The sym- 
plectic formulation of the problem allows us to write the per- 
turbation of the distribution function in terms of the gener- 
ating function g: this is equation ([ij. From this relation one 
can obtain the density 



p(q) 



PO+P 



(1) 



m/dp 

'm-{g,fS}d^P 



m/o d p ^ 



m/o d^p + 5pE + 5pL2 



For sufficiently small values of a, SpE is negligible in 
front of 5pj^2, as we have seen in ((20)) . hence (j'^'p = SpL2. 
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Using (|19|) . it can be checked that the first-order variation of 
totai mass associated to the perturbation is vanishing. 



1(1) 

m = 



5'^'pd'^q = m / 5pi^2 r'' s\n(6)Ar AO ^ 



dg 
dpe 



dpr dr dO c 



= mJ^iE) (^sinwJ^+cosW 
= 

Without more hypotheses than H23p on the perturbation 
generating function g, and for sufficiently small values of 
a, the first order induced variations of density are 



x(i) X 



f d^g ^ cosje) dg 
\d9dpe sin(6) dpe 



dpr 



In order to obtain some physical characteristics of the 
instability, we have to make yet another assumption about 
g. To find a function that verifies condition (|23p . given the 
form of g given in (|12p . we can suppose for example that g 
is separated in the 6 variable, i.e. one can find two functions 
A and B such that 

g{E,L\e,pe) = B{E,L\pe)A{e) (25) 

Under this assumption, criterion (|23|l becomes 



VS, 6» : A'(e)B(E,Q,Q) =Q while Ale)^lE,Q,Q) ^ Q 

ape 

(26) 

It is very easy to find a B that verifies this condition. A 
direct calculation then gives 



0^ ' p = m 



D(e) 



<^(i5)L2=o|^(^,0,0)dp. 



where 



If D 



is not constant, which corresponds to a wide 
class of then 5^^^ p does depend on and the spherical 
symmetry of the equilibrium state is broken.0 

We have reached our goal: we have found a class of 
perturbations g which leads to a negative energy variation, 
and which creates a density variation that is not spherically 
symmetric. As per section[2l this means that with the help of 
dissipation, the system is unstable against this perturbation: 
hence a favoured direction will appear in the system, which 
was initially spherical. 

Result 2. Consider a self-gravitating system, described by 
the CBP system and represented by a distribution function 
}o{E, L^), that is spherically symmetric and with nearly ra- 
dial orbits. Assume this system can dissipate energy. 

^ The equation D {8) ^ k can be easily solved and gives 

A(e) + h^l^ 

sm(0) 

wliere k and A are 9-free constants. 

^ The fact that the resulting perturbation depends only on 6, r 
and E, and thus is axisymmetric around the z axis, is of course 
a consequence of our choice of the form l|12p for the generating 
function. 



Then there exists perturbations, generated by a func- 
tion g, against which the system is unstable, and that cause 
it to lose its spherical symmetry. 



' 4 CONCLUSION 

In this paper we have shown two important points: self- 
gravitating dynamical systems described by the Colli- 
sionless Boltzmann-Poisson equations are candidates for 
Dissipation-Induced Instability when they are more and 
more radially anisotropic; and this mechanism generically 
introduces a favoured direction in the spatial part of the 
system's phase space. In comparison with previous tedious 
normal modes techniques used in this context, the detail 
of the first point gives a simple proof of radial orbit insta- 
bility based on energetics arguments. Dissipation, which is 
needed in our proof, is also implicitly required in the clas- 
sical intuitive understanding of this instability presented in 
iPalmerl l|l994 ) (section 7.3.1). It is a fact that in a pure ra- 
dial system — which is the most unstable — orbits, which 
are frozen in a fixed direction, cannot precess or librate as 
it is required for the trapping resonance invoked by Palmer. 
Hence, if two radial orbits actually attain a lower energy 
state by approaching each other, this mechanism actually 
needs a way to dissipate excessive energy. Finally, a point 
about time-scales should be stressed: it is well known that 
radial orbit instability is effective on a few crossing times, 
therefore if dissipation appears to be the cornerstone of ra- 
dial orbit instability, it is clear that it could not act alone. 
Non-linear and non-local aspects of the gravitational poten- 
tial clearly amplifies and completes the dissipation-triggered 
work. 
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